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Abstract. Models of self-organized criticality which can be described as sin- 
gular diffusions with or without (multiplicative) Wiener forcing term (as e.g. 
the Bak/Tang/Wiesenfeld- and Zhang-models) are analyzed. Existence and 
uniqueness of nonnegative strong solutions are proved. Previously numerically 
predicted transition to the critical state in 1-D is confirmed by a rigorous proof 
that this indeed happens in finite time with high probability. 

1. Introduction 

Within the past twenty years the notion of self-organized criticality (SOC) has 
become a new paradigm for the explanation of a huge variety of phenomena in 
nature and social sciences. Its origin lies in the attempt to explain the widespread 
appearance of power-law-statistics for characteristic events. In this paradigm an 
external perturbation may induce a chain reaction or avalanche in the system. Fur- 
thermore, a stationary state, the SOC-state, is reached where the average incoming 
flux is balanced by the average outgoing flux. This phenomenon was quite unex- 
pected since attaining the critical state of a thermodynamic system usually requires 
a fine tuning of some control parameter, which is absent in the definition of the 
SOC models. Several models have been proposed to mimic this mechanism includ- 
ing the sand pile BTW-model [TJ [2] and the Zhang- model [22] . The presence of 
thresholds in the definition of the dynamics implies that the energy can be accu- 
mulated locally, eventually generating a chain reaction which may transport energy 
on arbitrary large scales. 

The literature on SOC is vast. We refer e.g. to [TJ|2], [3], [7], [15], [3J, [PJHU], 
OH, DH, [13]: H]> Q3]> US, US], [22] for various studies. In [3J it was proposed to 
describe this phenomenon, e.g. in the case of the avalanche dynamics in the BTW- 
(see [HE]) and Zhang- (see [22]) models, by a singular diffusion. In the absense of 
noise the density g(t, £), t > 0, £ S R d , of this diffusion is formally described by the 
evolution equation 

(1) J^(*,0 = A¥(e(t,e)), 

where ^f(g) ■= f(g)H(g — g c ), g c > 0, H is the Heaviside function and f(g) = const, 
in the BTW-model and f(g) = g in the Zhang-model. Setting D(g) = f'(g)H(g — 
Qc) + f(g)S(g — g c ), the equation turns into 

(2) J^(t,0 = v-[P(fi(t,O)Ve(i,0] 

= D'(g(t, 0) |V<?(t, Of + D(g(t, £)) Ag(t, 0. 
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To discuss the problem in a heuristic way, let us consider a smooth version of H, 
for example 

H £ (g) = \ + iarctan(f) 

(or the mathematical more convenient one in ([7|) below) and the corresponding 
function D £ = D £ (g). Since H £ (g — g c ) is convex left of g c and concave right of 
Q c , we see that e.g. in the BTW-model D £ (g)\Vg\ 2 is positive if g is below g c and 
negative if g is above g c , i.e. according to g is "pushed" towards the critical 
value g c . This has been predicted numerically in 1-D by Bantay and Janosi in [3]. 

In [13] (see also [14]) Diaz-Guilera pointed out that it is more realistic to consider 
(Jl} perturbed by (an additive) noise to model a random amount of energy put into 
the system varying all over the underlying domain. So, ([I]) turns into a stochastic 
partial differential equation (SPDE). In [TT] (see also the references therein) based 
on numerical tests, Carlson and Swindle observed that in the presence of such a 
noise the self-organized behaviour does not necesssarily occur, i.e. the system fails 
to converge to the critical value g c . 

The purpose of this note is to provide rigorous mathematical proofs for the above 
phenomena. First, we sketch the proofs for existence and uniqueness of solutions to 
([!]) perturbed by noise (more precisely for multiplicative noise, so that positivity of 
initial data is preserved) . Second, we prove that at least in 1-D we have convergence 
to g c in finite time in the deterministic case (confirming the numerical results in 
[3]) and convergence to g c with high probability in the stochastic case. In regard 
to [11] one can probably not achieve more, but so far we failed to prove that this 
probability is really not equal to 1. 

Let us introduce our framework, where we switch to common notation in SPDE, 
i.e. replace g(i,£;) by X(t 7 x)(£) with x being the density at t = 0, t > 0, £ G O, 
where O is an open bounded domain in R d , d — 1,2, 3, with smooth boundary dO. 
The appropriate class of SPDE is then of the form 

!dX(t) - AV(X(t))dt 3 a(X(t))dW(t), 
in (0,oo) x O, 
*(I(i))9 0, on (0, oo) x dO, 
X(0,x) = x onO, 

where x is an initial datum, \P : M — > 2 R a maximal monotone graph and 

N 

<r(X)dW = VkXdfa e k , t > 0, 

k=l 

is a random forcing term, where {ek} C L 2 (0) is the eigen basis of the Laplacian 
—A on O with Dirichlet boundary conditions, N £ N U {+oo}, /i^ are positive 
numbers and (3k independent standard Brownian motions on a filtered probability 
space (O, J 7 , {f f }t>o,P). 

Throughout this note we make the following assumptions: 

Hypothesis 1. (i) *(r) = p sign r + *(r), for r G K, where G *(0) , p > 0, 
^ : M — > K is Lipschitzian, $£C (1\ {0}) and for some S > it satisfies 
$'(r) >S for allr GK\{0}. 
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(ii) If N = oo, the sequence {^k} is such that 



oo 



where are the eigenvalues of —A. 



A typical example is given by 

*(r) = tp (r) +c, r£l, 



where 



!a\r, r > 

r = 
a^r, r < 



and ai, a2 > 0, g > 0, c G \—g, g] arc constants. 

The following notations will be used. L p (0), p > 1, is the usual space of p- 
integrable functions on O with norm | • \ p . The scalar product in L 2 (0) and the 
duality induced by the pivot space L 2 (0) will be denoted by (•, -) 2 . Hq{0) C L 2 (0) 
is the first order Sobolev space on O with zero trace on the boundary. For a fixed 
measure space (E,£,m), a Banach space B and p £ [l,oo] we denote the space of 
all with respect to the measure m p-integrable maps from E to B by L P (E;B). 

In the following by H we shall denote the distribution space 



Jo 

We recall that the operator x — > — A$f(x) with the domain 

{x G ^(O) fl H : 3n e H^{0) s.th. r/ G ty(x) a.e. on O}, 

is maximal monotone in H (see e.g. [4]) and so the distribution space H is the 
natural functional setting for equation ([3]). However, the general existence theory 
of infinite dimensional stochastic equations in Hilbert space with nonlinear maximal 
monotone operators (see e.g. [12], [20]) is not directly applicable and so a direct 
approach must be used. 



Definition 2. Let x G H . An H -valued continuous Tt-adapted process X = X(t, x) 
is called a solution to ([3]) on [0,T] if for some p G [l,oo[ 





2. Existence, uniqueness and positivity 



X G L p (il x (0, T) x O) n £ 2 (0, T; L 2 (fi, iJ)), 



and there exists r) G L p (il x (0,T) x 0) such that P-a.s. 
(X(t,x),e 3 ) 2 

m = («iej}2 + / t /cj^(s J Ae i(0 c ^ ds 

V j G N, * G [0, T 



(5) 



r? G a.e. in x (0, T) x C 
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Below for simplicity we often write X(t) instead of X(t, x). 

Theorem 3. Under Hypothesis Q] for each x G L p (0), p > 4 there is a unique 
solution X to ([3]). Moreover, if x is nonnegative a.e. in O then P-a.s. 

X(t, x)(£) > 0, for a.e. (t, £) G (0, oo) x O. 

Sketch of Proof. Consider the approximating equation 

f dX x (t)-A(^ x {X x {t))dt = a(X x (t))dW(t), 
( ' 1 X x (0,x)=x, 

where A > 0, 

(7) $ A (r) := p (sign) A (r) + tf(r), rel, 

( 1 if r > A 
(sign) A (r) := J J ifr€[-A,A] 
[ 1 if r < -A. 

By [SI Theorem 2.2] (applied with m = 1) equation (J5]) has a unique solution 

X x G L 2 (ft x (0,T) x 0) nL 2 (O,C([0,T];iJ)) 

in the sense of Definition [2] which is nonnegative, if so is x. Here as usual the space 
of continuous H- valued paths C([0, T];H) is equipped with the supremum norm. 

By Ito's formula for a > large enough it follows that for all A, fj, G (0, 1) and 
t G [0, T] 

\\X x {t)-X,{t))\ 2 -ie-«t 
<Cmax{A,M}/o/e,(l*A(X A ( S ))| 2 + |X A ( S )| 2 

+ |^(^( S ))| 2 + |X M ( S )| 2 ) e —^^ 

+ f*e-<* s {X x (s) - X^(s),a{X x (s) - X f ,( s ))dW( S )) 2 . 

Hence by the Burkholder-Davis-Gundy inequality {X\} is a Cauchy net in L 2 (Q; C([0, T], if)) 
and by a standard technique from stochatic PDE one shows that the limit X is the 
desired solution to (cf. [S]). □ 

Remark 4. One can also show (see Prop. 3.4]; thatX, X x G L 2 (Q, T: L 2 (il, Hi (£>))), 
i/iaf 



limfi/ |X A -X| 2 2(o) dt = 0, 



where E denotes expectation with respect to P, and £/ia£ &o£/i X and X A /lave con- 
tinuous paths in L 2 (0). Theorem^ is true for more general not linear growing ^> 
(see W- 

3. Extinction in finite time and self-organized criticality 
In this section we assume N < oo. Let r be the stopping time 
r = inf{t > : \X(t,x)\-i = 0}, 
where X(t, x), t > 0, is the solution from Theorem[31 
Proposition 5. 

X(t,x) = for t > t, P-a.s.. 
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Sketch of Proof. For simplicity we consider the case with g = 1. Define 

N N 

Kt) '■= -^VkCkPkit), te[0,T], fi :=^2n 2 k el 

k=l k=l 

and 

Y(t) := e^ (t) X(i), t > 0. 
Then by Ito's product rule Y satisfies P-a.s. the following ordinary PDE 

^=e^Ar,(t)-^Y(t), t>0. 
Setting Y\ :— e^X x we consider the approximating equation 

(9) ^) =e Kt) A ^ (t) _l /i(t)n(t)) f > 0j 

where 

(10) 7 7A (t) = * A (X A (i))Gif 1 (O). 



Hence 



(ii) 



dt 

= (r lx (t),A(e^Y x (t))) 2 -^(mYx(t),Y x (t)) 2 , 
where by (fTP)) , and integrating by parts we have 



A 



( m (t),A(e^Y x (t))) 2 

' / (IvnWI 2 

/{e-M*)|y A (t)|<A} 

-|y A (t)| 2 |W)I 2 R 
§'(e-"Wy A (t)) 

(|vr A (i)| 2 -|y A (<)| 2 |V/i(<)| 2 )^. 

This yields 

(12) ( m (t),A(e^Y x (t))) 2 <C(\Y x (t)\l + X). 
Hence (If 1 [) and Gronwall's lemma imply 

Now letting A — > we get 

(13) \Y(t)\j <e c ^\Y(s)\l t>s. 

Taking in (|f 3|) s = twc get Y"(i) — for all t > t as claimed. □ 

For proving the extinction result we need 0CM, i.e. d = 1. To be more specific 
let = (0,?r). Then e fc (£) = sinfcf, £ £ [0,tt], A fc = fc 2 and 1^(0, tt) C H 
continuously, so 



(14) 7 



inf/M^ : x G (0,tt)1 > 0. 
L|x|_i J 
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Theorem 6. Consider the equation 

' dX(t) - A(p sign (X(t) - x c ) + V(X(t) - x c ))dt 

N 

3 a(X(t) - x c )^2fi k e k d/3 k , t>0, 



(15) 



fe=i 

p sign (X(t) - x c ) + <&(X(t) - x c ) 3 0, on 9(0, n), 
. X(0,x) = x. 



where i c eR. 
Assume that 

\x - x c \-i < pjC^ 1 , 
where Cn := \ SftLi(l + ^) 2 / i fe an d 7 * s as * n GM- Then for all t > 

(16) P(r,<t)>l- ''" g ^ , 

t c = inf{t > : \X(t) - x c |_i = 0} 
= sup{t > : - x c |_i > 0} 

and X — x c is the solution from Theorem^ 

Sketch of Proof. For simplicity we assume x c = 0. Since (X (t) , ^ (X (t))) 2 > 0, 
6 *(0). An application of Ito's formula for ip^XW) = {\X\ 2 _^ + e 2 ) 1/2 yields 
P-a.s. 



min(£Ar 



1*001- 



(|A»| 2 _ 1+£ 2)l/2' 

2 I (a(X(s))dW(s),^(\X(s)\ 2 _ 1 )X( S )). 



Now, letting e tend to zero we get 'P-a.s. for t > 

+ 7pmin(i A r) 

( 18 ) < M-i +Cjv/o \X(s)\-ids 

+2 / l [0iT] ( S )(cr(X( S ))dW( S ),X( S )|X( S )|:J). 
Jo 

Hence by a standard comparison result 

|*(*)|-i+7P / e CN( ^ s) l[o,r](s)ds 



<e c "*|x|_ 1 + / HX( S ))dW^s),A%s)|X( S )|:J>. 
Jo 

Taking expectation we get 

V CnS P(t > s)ds < 

7P 

Writing P(r > s) = 1 - P(r < s) we deduce (US). □ 
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Corollary 7. If in the situation of the above Theorem^ the noise is zero, i.e. 
Cn = 0, then 

\x- X c \-x 

T c < • 

Ql 

Proof. The assertion follows from (|18p . □ 
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